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ON  11  CM  IRREGULARITIES  DURING  EQUATORIAL  SPREAD  F 
I.  Introduction 

During  the  past  several  years,  high-frequency  radar  backscatter 
experiments  have  revealed  a  spectrum  of  short -wavelength  (i.e.,  below  the  ion 
gyroradius)  irregularities  during  equatorial  spread  F  (ESF) .  Radar  back¬ 
scatter  observations  at  50  MHz,  155  MHz  and  415  MHz  indicate  density 
fluctuations  exist  with  scale  sizes  of  3m,  lm  and  36  cm,  respectively  [Farley 
et  al . .  1970;  Woodman  and  LaHoz,  1976;  Costa  and  Kelley,  1978a, b;  Huba  et  al . , 
1978].  Most  recently,  Tsunoda  (1980)  has  observed  radar  backscatter  from 
11  cm  (1320  MHz)  irregularities  during  equatorial  spread  F  at  high  altitudes, 
using  the  TRADEX  radar.  These  observations  were  part  of  a  coordinated  Defense 
Nuclear  Agency  campaign  at  Kwajalein  to  study  ionospheric  irregularities 
during  equatorial  spread  F.  The  program  was  designed  to  obtain  simultaneous 
radar  backscatter  results  and  in  situ  rocket  measurements  of  density  and 
fluctuating  fields.  Unfortunately,  the  11  cm  observations  were  made  prior  to 
the  rocket  launch  (but  during  ESF) ,  so  coincidental  in  situ  plasma  and 
fluctuation  data  are  not  available.  However,  presumably  ionospheric  plasma 
conditions  were  similar  during  the  various  measurements. 

Sharp  density  gradients  were  observed  during  this  campaign  (M.  C. 

Kelley,  private  communication,  1980)  and  have  been  observed  during  past 
equatorial  spread  F  events  (Costa  and  Kelley,  1978a, b).  The  scale  lengths  of 
these  gradients  range  from  tens  of  meters  to  several  hundred  meters  and  are 
presumably  due  to  primary  longer  wavelength  instabilities  such  as  the 
Rayleigh-Taylor  instability.  Since  the  typical  ion  gyroradius  is  rLi  ~  5m, 

it  is  found  that  rT./L  ^0.2  where  L  is  the  density  gradient  scale  length. 

Li  n  n 

Based  upon  this  evidence,  it  has  been  suggested  that  various  drift  instabi¬ 
lities  are  responsible  for  the  short  wavelength  irregularities  [Huba  et  al. , 
1978;  Costa  and  Kelley,  1978a, b;  Huba  and  Ossakow,  1979a, b],  depending  upon 
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the  wavelength  observed.  Although  collisionless  drift  waves  would  easily  be 
excited  under  these  circumstances,  collisional  effects  have  a  stabilizing 
Influence  on  the  instabilities  investigated  thus  far  [Huba  and  Ossakow. 

1979a, b;  Sperling  and  Goldman.  1980],  Specifically,  the  lower-hybrid-drift 
instability  is  the  prime  candidate  to  explain  the  lm,  36  cm  and  11  cm  irregu¬ 
larities.  However,  recent  work  [Huba  and  Ossakow.  1979a;  Sperling  and 
Goldman.  1980]  has  indicated  that  electron-electron,  electron-ion  and 
electron-neutral  collisions  stabilize  the  instability  and  an  approximate 
threshold  condition  has  been  derived  [Sperling  and  Goldman.  1980]. 

In  the  light  of  the  new  observations  of  11  cm  irregularities  and 
theoretical  results,  we  show  that  the  lower-hybrid-drift  instability  is  the 
most  probable  cause  of  the  small-scale  irregularities  (i.e.,  11  cm,  36  cm  and 
lm)  observed  during  equatorial  spread  F.  Moreover,  we  discuss  the  role  of 
these  irregularities  on  the  large  scale  plasma  behavior  during  these  periods. 

II.  Theory 

We  consider  a  plasma  immersed  in  a  homogeneous,  unidirectional  field 

JB  *  Bq  ez  with  an  inhomogeneous  density  profile  n  ■  n(x)  and  a  constant 

temperature  (Tg  ■  *  T).  Each  species  a  (i.e.,  electrons  and  0+  ions)  has 

a  diamagnetic  drift  velocity  ■  (va2/2fia)  d  In  n/dx  where  va  -  ^Tg/n^)1* 

is  the  thermal  velocity,  Qa  ■  eg^Q^mo,c  is  the  cyclotron  frequency  and 

n  «  n  «  n . .  A  net  current  exists  in  the  plasma  J  ■  en(V . .  -  V.  )e 
e  i  dl  de  y 

2en  ey  which  provides  the  free  energy  to  drive  an  instability.  Perturbed 
quantities  are  assumed  to  vary  as  exp[iky  -  luut]  and  we  consider  only  electro¬ 
static  oscillations  since  [S  «  1.  We  make  use  of  the  local  approximation 

which  requires  kL  »  1  where  L  ■  (d  in  n/dx)-1  is  the  density  gradient 
n  n 

scale  length.  The  ions  are  taken  to  be  unmagnetized  while  the  electrons  are 
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considered  to  be  magnetized.  The  ions  behave  as  unmagnetized  particles  for 
the  wavelengths  under  consideration  (k  r^  ~  102)  because  of  ion-ion  col¬ 
lisions  (Huba  and  Ossakow.  1979a).  Electron-neutral,  electron-electron  and 
electron-ion  collisions  are  included  in  the  analysis  via  an  effective 
collision  frequency  vg  -  vgn  +  vgl  (1  +  0.15  k2rLfi2)  where  r^  ■  ve/|Qel  is 
the  mean  electron  Larmor  radius  (Kadomtsev.  1965;  Sperling  and  Goldman,  1980). 
We  note  that  our  coefficient  of  k2r^g2  differs  from  that  of  Sperling  and 
Goldman  (1980)  due  to  a  different  definition  of  the  thermal  velocity. 

Based  upon  these  assumptions,  we  obtain  the  following  dispersion 
equation  for  the  lower-hybrid-drift  Instability 


D(tu,k)  -  1  +  xt  +  Xfi  (1) 

where  (Huba  et  al. ,  1978) 


where  V2  “  4™aea2/V  va2  *  2VV  na  "  WV*  vda  *  (va2/2CVd  **  n/dx* 
Ve  “  ven  +  vei(1  +0-3be)’  be  *  ^u^2*  o(x)  "  1n  is  the  ""dified 

Bessel  function  of  order  n  and  Z  is  the  plasma  dispersion  function.  We  note 

that  Eq.  (3)  is  based  upon  the  BGK  collision  model  and,  strictly  speaking, 
does  not  correctly  treat  electron-electron  collisions  (Rukhadze  and  Sllln. 
1968).  However,  electron  viscosity  is  approximately  modeled  via  the  term 
proportional  to  k2r^e2  (Mikhallovskli  and  Poautse.  1966).  Since  the  11  cm 
irregularities  correspond  to  kr^e  *  2 ,  it  is  clear  that  electron  vlacotw 
effects  are  only  moderately  Important.  Moreover,  the  BGK  model  is  adequate 


in  the  absence  of  temperature  gradients  (Rukhadze  and  Silin.  1968)  which  is 
the  situation  in  the  F  region.  Thus,  Eq.  (3)  can  be  used  with  confidence  to 
describe  the  electron  response  qualitatively.  The  quantitative  results  based 
upon  Eq.  (3)  are  approximately  correct  since  a  model  Fokker  Planck  equation 
is  used  to  describe  the  collisionality . 


Since  «  v^  (which  corresponds  to  r^  «  Ln),  we  can  expand  the 
plasma  dispersion  function  in  the  small  argument  limit  (i.e.,  Z(t| »)  i/iT). 
The  dispersion  equation  assumes  the  form 


D(ui,k) 


^pi2  r 

1  +  k^  L 


i  +  isa 


2w  2 
pe 


k2v  2 
e 


uu-kV .  +  iv 

de  e 


r  (b  )  i 


w] 


In  the  limit  y  ~  vg  «  uj^.,  where  u>  -  uur  +  iy  ,  Eq .  (4)  has  the  following 


solution 


•r  ■  ro  kvdi  [*  +  k2>di2+  <VV(l-rj 


■/-  UJ  -kV . .  T.  V  /  kV  T  \' 

-  -  Wk'V  [r  -  w7~ +  r  jr  <l-r„>(1+  -if  f) 


where  X^2  “  v^2/2u)pjf  and  the  argument  of  has  been  suppressed.  In  the 

absence  of  collisions  (v  -  0),  instability  occurs  for  w  <  kV..  and  there  is 

e  r  di 

no  threshold  requirement.  However,  electron  collisions  are  stabilizing  and 
place  a  threshold  condition  on  the  drift  velocity  to  excite  the  mode  (Huba 
and  Ossakow,  1979a; Sperling  and  Goldman,  1980).  Substituting  Eq.  (5)  into 
Eq.  (6),  the  critical  drift  velocity  (i.e.,  such  that  y  >  0)  is  given  by 


/vdi\  (i-v  -  I0)2<2+k2V>T! 

V«i /  Mv  /r  k^„  +  k2>di2  . 


4 


This  corresponds  Co  a  critical  density  gradient  scale  length  via 

<  <VVdi>cr  <rL, 

We  now  apply  Eq.  (7)  to  plasma  conditions  relevant  to  equatorial 
spread  F  to  determine  whether  or  not  the  lower-hybrid-drift  instability  is 
responsible  for  the  11  cm  radar  backscatter  observations.  We  choose  B  =  0.3  G, 
Tg  *  *  0.1  eV  and  m^  ■  16  n^.  The  collision  frequency  is  given  by 

v  -  v  +  v  .(1  +  0.15  k^r.  2)  (8) 

e  en  ei'  Le  '  '  ' 

where  (Braginskii.  1965;  Johnson,  1961) 

vm  *  5.0  x  10-8  sec-1  (9) 

vel  «  (X/3.5  x  105)(ne/Te3/2)  sec-1  (10) 

and  nn  is  the  neutral  density,  ng  is  the  electron  density 

X  =  23.4  -  1.15  log  n^  +  3.45  log  Te,  and  Tg  is  given  in  eV.  We  plot 

(Vdi/vi)  and  Lncr(meters)  versus  krLg  in  Fig.  1  for  the  following  values: 
cr 

*  104,  105,  106  cm  3  and  nn  «  10®,  5  x  10®  cm  3.  The  range  in  electron 

density  corresponds  to  the  ambient  density  (ng  ~  106cm-3  )  to  the  density 

within  ionospheric  bubbles  or  plasma  depletions  (ng  ~  104cm  3).  The  values 

of  the  neutral  density  are  for  altitudes  of  500  km  (n  ~  10®cm  3)  and  240  km 

n 

(nR  ~  5  x  10  cm-3)  during  sunspot  maximum  (Johnson,  1961).  The  altitude  of 
500  km  corresponds  to  the  region  from  which  the  11  cm  radar  backscatter 
returns  were  observed  (Tsunoda,  1980);  whereas,  the  240  km  altitude  corres¬ 
ponds  to  the  bottom  of  the  F  region  on  that  same  night.  Also,  we  have 

marked  the  values  of  kr,  which  correspond  to  lm,  36  cm  and  11  cm,  res- 

Le 

pectively.  Several  interesting  features  of  Fig.  1  are  as  follows. 
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First,  minimum  values  of  (V,./v. )  and  L  cr  exist  for  kr.  ~  0.6  as 

ai  l  cr  n  Le 

noted  by  Sperling  and  Goldman  (1980) .  The  minimum  is  rather  broad  for  low  n 
-  e 

but  becomes  sharper  as  nfi  is  increased.  This  indicates  that  for  a  given 

density  and  density  gradient  scale  length  only  a  certain  range  of  kr,  can  be 

Le 

excited  linearly.  Second,  larger  values  of  the  neutral  density  require 

larger  drift  velocities  (or  shorter  density  gradient  scale  lengths)  as 

expected.  This  is  a  more  dramatic  effect  at  lower  electron  densities 

(ng  ~  104cm-3)  than  higher  electron  densities  (ng  ~  106cm“3)  simply  because 

v  ,/v  scales  as  n  .  And  finally,  for  n  ■  10**cm“3  there  is  a  cutoff  at 
ei  n  e  e 

kr.  1.5.  This  arises  because  kX.  »  1  for  these  values  and  one  does  not 
Le  de 

expect  collective  plasma  phenomena  to  occur  on  length  scales  shorter  than  a 
Debye  length.  Thus,  for  electron  densities  less  than  104cm“3,  there  will  be 
no  density  fluctuations  with  scale  lengths  shorter  than  20  cm  other  than 
thermal  fluctuations. 

III.  Discussion 

Recently,  during  a  coordinated  ground-based  rocket  campaign  to  study 
ionospheric  irregularities  during  equatorial  spread  F  at  Kwajalein,  Tsunoda 
(1980)  observed  radar  backscatter  from  11  cm  irregularities.  These  are  the 
smallest  scale  irregularities  observed  thus  far  and  are  comparable  to  the 
mean  electron  Larmor  radius  (rLe  ~  3  cm) .  An  important  feature  of  the 
experimental  results  is  that  these  irregularities  were  only  observed  at  high 
altitudes  (~  500  km).  During  the  rocket  flight,  which  occurred  subsequent  to 
the  II  cm  radar  backscatter  measurements,  in  situ  probes  detected  sharp 
density  gradients  (Lr  £  45m;  Kelley,  private  communication,  1980)  and  large 
density  depletions  (ng  ~  several  x  104cnf3;  E.  Szuszczewicz,  private  communi¬ 
cation,  1980)  at  high  altitude  (~  500  km).  We  note  that  previous  in  situ 
rocket  measurements  have  found  density  gradient  scale  lengths  as  small  as  30m 
(Costa  and  Kelley,  1978a). 
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These  data  suggest  the  following  scenario.  As  equatorial  spread  F 

develops,  density  depletions  rise  to  the  topside  of  the  F  region  where  the 

neutral  density  is  low  (~  500  km).  Within  these  plasma  bubbles  (where 

electron  collisional  effects  are  minimal)  sharp  density  gradients  exist  which 

can  excite  the  lower-hybrid  drift  instability.  The  density  fluctuations 

associated  with  the  instability  give  rise  to  the  radar  backscatter  measure- 

5  -3 

ments  at  lm,  36  cm  and  11  cm.  This  can  occur  for  L  20m,  n  <  10  cm  and 

n  *  e 

8  ~3 

n  ~  10  cm  (see  Fig.  1)  which  is  consistent  with  observational  results  to 
n 

within  a  factor  of  two.  This  is  quite  good  considering  the  approximations 

involved  in  deriving  the  electron  response  (Eq.  (3)).  At  lower  altitudes 

(<  250  km)  the  increased  neutral  density  requires  a  sharper  density  gradient 

scale  length  for  a  given  value  of  n  (or  lower  values  of  n  for  fixed 

e  e 

density  gradient  scale  length)  especially  for  the  11  cm  irregularities 

(see  Fig.  1)  and  the  mode  is  stable.  Also,  it  is  possible  that  the  density 

4  -3 

within  the  lower  altitude  density  depletions  is  ng  10  cm  and  radar  back¬ 
scatter  is  not  observed  at  11  cm  since  this  is  smaller  than  the  Debye  length. 
Thus,  based  upon  the  linear  theory  of  the  lower-hybrid-drift  instability,  it 
is  found  that  this  mode  is  the  most  probable  cause  of  the  small-scale 
irregularities  (^  lm)  observed  during  equatorial  spread  F. 

Finally,  we  discuss  the  role  of  this  instability  in  the  evolution  of 
the  large  scale  plasma  phenomena  occurring  during  equatorial  spread  F. 
Typically,  plasma  microturbulence  influences  the  plasma  via  its  anomalous 
transport  properties  (i.e.,  scattering  of  particles  by  the  collective  fields 
associated  with  the  instability).  A  recent  theoretical  study  of  the  lower- 
hybrid-drift  instability  (Drake ,  1980)  Indicates  that  this  mode  can  only 
produce  irreversible  electron  heating  and  diffusion  when  e$/T  ^  0.2-0. 5 
where  $  is  the  fluctuating  electrostatic  potential.  To  produce  this  level 
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of  turbulence  requires  very  sharp  density-gradients  (L  ^  r  ,  )  which  do  not 

n  JjX 

exist  during  equatorial  spread  F.  For  typical  ionospheric  conditions  (i.e., 

L  »  rT .),  the  Instability  will  saturate  at  a  low  level  of  turbulence  and 
n  Li 

not  be  an  effective  anomalous  transport  mechanism.  Thus,  the  small-scale 
density  and  field  fluctuations  associated  with  this  mode  are  a  signature  of 
equatorial  spread  F  (under  the  proper  conditions)  and  will  probably  not  sig¬ 
nificantly  influence  the  macroscopic  fluid  evolution  of  the  plasma  for  scale 
sizes  Ln. 
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1m  36cm  krLe  11cm 

Fig.  1  —  Plot  of  (Vdi/Vi)cr  and  Iff  (meters)  vs.  k  r^.  Values  of 
Vdi/v j  and  L„  above  the  curves  lead  to  instability,  while  those 
below  lead  to  stability.  The  following  densities  are  considered: 
electron  density  (ne  =  104, 106, 106cm—  3)  and  neutral  density 
(n„  =  108cm-3  and  nn  =  6  X  109cm~3...). 
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